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When using distance measurements to probe spatial curvature, the geometric degeneracy between 
curvature and dark energy in the distance-redshift relation typically requires either making strong 
assumptions about the dark energy evolution or sacrificing precision in a more model-independent 
approach. Measurements of the redshift evolution of the linear growth of perturbations can break the 
geometric degeneracy, providing curvature constraints that are both precise and model-independent. 
Future supernova, CMB, and cluster data have the potential to measure the curvature with an accu- 
racy of a (Ok) = 0.002, without specifying a particular dark energy phenomenology. In combination 
with distance measurements, the evolution of the growth function at low redshifts provides the 
strongest curvature constraint if the high-redshift universe is well approximated as being purely 
matter dominated. However, in the presence of early dark energy or massive neutrinos, the preci- 
sion in curvature is reduced due to additional degeneracies, and precise normalization of the growth 
function relative to recombination is important for obtaining accurate constraints. Curvature limits 
from distances and growth compare favorably to other approaches to curvature estimation proposed 
in the literature, providing either greater accuracy or greater freedom from dark energy modeling 
assumptions, and are complementary due to the use of independent data sets. Model-independent 
estimates of curvature are critical for both testing inflation and obtaining unbiased constraints on 
dark energy parameters. 



I. INTRODUCTION 

Cosmological measurements of the average spatial cur- 
vature of the spacetime metric are one of only a handful 
of methods available for testing the inflationary paradigm 
for the early universe. Current observations are consis- 
tent with the inflationary prediction of a nearly flat uni- 
verse. However, the precision of curvature measurements 
is only at the percent level at best (e.g., [l|, 0), whereas 
the expected level of curvature in standard inflationary 
scenarios is generally much smaller. Moreover, obtaining 
the most precise limits on curvature requires assuming a 
particular simple form for the dark energy evolution due 
to the well-known "geometric degeneracy" between cur- 
vature and dark energyijl i, i, 0, 1, ©, El EH, El, El, 
E3, El, El, E3, El, El, M^M ■ Inferences about inflation 
based on such curvature constraints are only valid if the 
assumed dark energy behavior is an adequate description 
of the true evolution. 

Just as uncertainty about the dark energy evolution 
affects estimates of curvature, uncertainty about curva- 
ture limits our ability to constrain parameters of dark 
energy models with cosmic distances. One often assumes 
spatial flatness motivated by the predictions of inflation 
when constraining dark energy models, but the resulting 
parameter estimates may be biased if the true spatial cur- 
vature deviates even slightly from zero [1, EH, El, EH, EH • 

A few methods for using measured distances to obtain 
curvature estimates that are independent of the dark en- 
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ergy evolution have been proposed for use with future 
data sets. For example, Bernstein [l8[ proposed a tech- 
nique using weak lensing galaxy-shear correlations to 
measure triangles of distances between the lensing and 
source galaxy planes and the observer, resulting in a cur- 
vature estimate that depends only on the assumed form 
of the spacetime metric. An alternate method by Knox 
[l9t uses precisely measured distances at high redshift 
( z si 3) combined with the distance to recombination 
from cosmic microwave background (CMB) data to in- 
fer the curvature without dependence on the low-redshift 
dark energy evolution. 

In this work, we describe a different approach to model- 
independent curvature estimates that uses combinations 
of data sets that probe the distance-redshift relation 
and the growth of linear perturbations. The geomet- 
ric degeneracy in distance data arises because distances 
depend on both the expansion rate (and therefore the 
dark energy evolution) and the spatial curvature. The 
growth of structure, on the other hand, depends only on 
the expansion rate. If one can measure both distances 
and growth over a similar range of redshifts, then con- 
straints on the expansion rate from growth data can be 
used to break the degeneracy in distance data, provid- 
ing a model-independent determination of the curvature. 
The only assumptions required are that general relativ- 
ity (GR) is the correct theory of gravity governing the 
growth of structure and that dark energy does not clus- 
ter significantly on the scales of interest. 

In Section [HI we review the basic distance and growth 
relations and observables. Section [TTTl describes the geo- 
metric degeneracy in distances and how this degeneracy 
is broken by growth information. We then present fore- 
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casts for curvature constraints from future distance and 
growth data, beginning with descriptions of two meth- 
ods of relating growth evolution on linear scales to an 
observed distance-redshift relation. The first, in Sec. lIVl 
is based on a numerical exploration of general dark en- 
ergy models carried out in Ref. [12] by Mortonson, Hu, 
and Huterer (hereafter, MHH09). The second method in- 
volves rewriting the equations for distances and growth 
so that the common dependence on the expansion rate 
drops out. The basis of this method of reconstructing 
the growth history from observed distances comes from 
Alam, Sahni, and Starobinsky [23|, and in Sec. [V] we 
summarize this work and extend it in several ways to 
allow the method to be applied to curvature forecasts. 
These two methods are complementary in several ways; 
the MCMC approach is more straightforward in terms of 
error propagation and the solution for the growth evo- 
lution, but it can be quite time-consuming and depends 
more on one's priors on the dark energy evolution com- 
pared with the analytic growth reconstruction approach. 
The growth reconstruction method is therefore a useful 
tool for exploring the curvature-dependent relation be- 
tween distances and growth for a variety of different cos- 
mologies and assumed data sets, while the MCMC re- 
sults help in testing and calibrating the analytic method 
and in providing accurate error estimates. Using both of 
these methods, forecasts for curvature from a combina- 
tion of future supernova (SN), CMB, and X-ray cluster 
data are presented in Sec. IVII Finally, Sec. IVIII contains 
a summary and discussion of the results of this work. 



II. PRELIMINARIES 
A. Spatial curvature 

Given that the universe appears to be spatially homo- 
geneous and isotropic on large scales, the background 
metric can be written in the Friedmann-Robertson- 
Walker (FRW) form: 



ds 2 



-dt z 



dD 2 



1 + n K H 2 D 2 



D 2 (d0 2 



(1) 

which describes an expanding (or contracting) universe 
with scale factor a(t), where D is the comoving radial 
coordinate and Hq is the Hubble constant. The FRW 
metric has constant spatial curvature parametrized by 
fix, where a flat universe has Qk = 0, an open universe 
f^K > 0, and a closed universe f^K < 0. The curvature 
parameter is related to the total density of the compo- 
nents of the universe in units of the critical density for 
flatness, tt to t = /O t ot/p C r,o, by fl K = 1 - O tot , where all 
densities are evaluated at the present time. 

Theories of inflation predict that the universe is nearly 
flat (f^K ~ 0), and the fact that current observations are 
consistent with flatness is viewed as supporting evidence 



for inflation. The expected deviations from flatness are 
typically at or below the level of the initial curvature 
perturbations at the end of inflation, |Ok| ^ 10~ 5 (e.g., 
[H HE])- The ultimate precision with which the cur- 
vature may be determined from observations is limited 
by cosmic variance and model selection considerations to 
a{n K ) ~ KT 5 - 1CT 4 [li,|23l. While there are some the- 
ories of inflation in which the present value of the curva- 
ture is large enough to be potentially observable without 
excessive fine tuning of the initial conditions of inflation 
HE HI HI El EE E3, El El El El El , a detection of 
nonzero curvature would challenge at least the simplest 
inflationary theories. 

The strongest observational bounds on curvature are 
presently at the percent level, i.e. ct(Ok) ~ 0.01. The 
main limits on curvature come from measurements of 
angular diameter distances in the CMB at z ~ 1000 
and BAO at z < 1. However, these constraints rely 
on assuming a simple model for the dark energy such 
as a cosmological constant. More precise and more 
model-independent measurements of the spatial curva- 
ture would provide valuable tests of theories of inflation. 



B. Distances 

In a flat universe, the comoving distance Df to an ob- 
ject at redshift z obtained by integrating over the comov- 
ing radial coordinate in the FRW metric is 



Df(z) 



dz' 

o wr 



(2) 



More generally, for universes with nonzero spatial curva- 
ture the comoving distance is 



D{z) = -S K [KD f (z)], 

K 



(3) 



where K = (IOkI-^o ) 1 ^ 2 ^ s the inverse of the curvature 
radius of the universe, and Sk(x) = x for a flat universe, 
sinhx for an open universe, and sin a; for a closed uni- 
verse. Distances therefore depend on both the expansion 
rate, H(z), and geometry, fix- The distances at low red- 
shifts for three models with varying spatial curvature are 
plotted in Fig. [1] Note that the curvature dependence is 
very weak at low rcdshifts, but high-redshift distances, 
e.g. the distance to recombination, are more sensitive to 
the geometry of the universe. 

The main probes of distances we will consider here are 
Type la supernovae (SNe), standardizable candles whose 
average magnitudes are related to distances as 



i(z) = 5log[H d h (z)] + M, 



(4) 



where M = M — 5 log(i/o/Mpc _1 ) + 25 combines the 
unknown absolute magnitude of the supernovae M and 
Hubble constant Hq, both of which only affect the overall 
normalization of the SN distance-redshift relation. Since 
it involves unknown parameters that do not affect the 
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FIG. 1: Comoving distance (top), growth relative to z — 
(middle), and growth relative to recombination (bottom) 
vs. redshift for flat, open, and closed models. For all three 
models, Q m — 0.24 and h = 0.73. In the top panel, the 3 
curves are indistinguishable in HqD(z). 



measured evolution of distances with redshift, M is a 
"nuisance" parameter that is generally marginalized in a 
cosmological analysis of SN data. 

Because the distance normalization is unknown, SN 
data determine relative distances, but not the abso- 
lute scale of the distance-redshift relation. However, 
SNe at low z can constrain the normalization since 
liin z _»rj HqD(z) = z under reasonable assumptions about 
the evolution of H(z) in the recent past. 1 Then for 
low-redshift SNe, the average magnitude from Eq. (H} 
is m(z) w 5 log z + M. which provides an estimate of A4. 

We will also consider the angular diameter distance 
constraint from the acoustic scale of the CMB. The main 
effects of dark energy and curvature on the CMB enter 
through the distance to recombination at z* as 1089 [H[ 
and the matter density fl m h 2 , which affect the angular 
scale and amplitude of the acoustic peaks [2(J. Dark 
energy and curvature also affect the large-scale CMB 
anisotropics through the integrated Sachs- Wolfe effect, 
but the information available is limited due to cosmic 
variance on those scales and the resulting constraints are 
relatively weak. We will therefore neglect this informa- 
tion in the curvature forecasts. 



1 For example, assuming that there was not a sudden large tran- 
sition in the dark energy equation of state at z < 0.01 (e.g., see 
IMHH090 



Angular diameter distances can also be determined 
by measuring baryon acoustic oscillations in the mat- 
ter power spectrum in the plane transverse to the line 
of sight. BAO can be a powerful probe of absolute dis- 
tances, but incomplete redshift coverage and the need 
for wide redshift bins make the technique less suitable 
than Type la SNe as a primary source of the distance in- 
formation for predictions of the growth evolution. How- 
ever, BAO data can provide complementary constraints 
on curvature through other means (e.g., [191]). 



C. Growth of linear perturbations 

The growth of linear matter perturbations obeys 

5 + 2H6-4wG N p ia 6 = (5) 

where S = Sp m / p m and overdots are derivatives with re- 
spect to time t. We assume here and throughout this 
work that general relativity is valid and that the dark 
energy is smooth on the relevant scales so that additional 
terms in the growth equation describing the clustering of 
dark energy can be neglected. 

Equation can be written in terms of G oc (1 + z)S 

as 



d 2 G 
din a 



din a 
dlnH 
din a 



dG 
din a 
3 



(6) 



G = 0, 



where fi m (z) = ri m H 2 (l + z) 3 /H 2 (z). The growth vari- 
able G is constant in a universe that contains only mat- 
ter, so it is nearly constant at high redshifts during mat- 
ter domination. We normalize the growth functions to 
S(z = 0) = 1 and G(z — > oo) = 1 so that 



(l + z)5(z) 



G(z) 



Go(z). 



(7) 



G(z = 0) 

Figure Q] shows Go(z) and G(z) for three models with 
different values of the spatial curvature. 

From Eqs. ((3]) and ((6]) one can see that distances de- 
pend on both the expansion rate and geometry, while 
growth depends only on the expansion rate (and Q, m H^ 
which is well determined by CMB data). Combinations 
of distance and growth information with similar redshift 
coverage therefore determine the geometry of the uni- 
verse with reduced dependence on the expansion rate. 

Measurements of cluster abundances in a range of red- 
shift bins can probe the growth evolution at low red- 
shifts, determining Gq(z). Clusters can also constrain 
G(z = 0) to probe high-rcdshift changes in the growth 
evolution due to massive neutrinos or early dark energy 
by comparing the growth determined by low-z clusters 
with the predicted growth extrapolated to low redshifts 
from measurements of the CMB power spectrum ampli- 
tude at z ~ 1000 (e.g., see [1^]). Ref. [5o] contains cur- 
rent examples of both types of measurements using an 
X-ray cluster sample. 
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Assuming that dark matter halos of mass M host 
galaxy clusters with the same mass, the cluster abun- 
dance depends on cosmology primarily through the halo 
mass function dn/dM and the comoving volume element 
in a solid element d£l and rcdshift slice dz, 



dV _ D 2 (z) 
dfl dz ~ H(z) ' 



(8) 



The mass function describing the comoving density of 
dark matter halos of mass M at redshift z can be written 



as 



dm 
dM 



^mPcr,o din a 1 (M,z) 



M 



dM 



f(a(M,z)). (9) 



Here o 2 (M, z) is the variance of linear matter density 
perturbations, 



cr 2 (M, z) = 



G (z) 



1 



dlnfc A 2 {k)W 2 (kR(M)), (10) 



where R(M) = [3Af/(47rfi m p CI , )] 1/3 , A 2 (fc) is the di- 
mensionless power spectrum of linear matter pertur- 
bations at z — with comoving wavenumber k, and 
W(kR) = 3ji(kR)/kR is the Fourier transform of a 
spherical top-hat window function with radius R. 

The function f(a) in Eq. §§§ parametrizes the mass 
function in a way that is relatively independent of red- 
shift and cosmological parameters. The dependence on a 
is exponential as a — > 0, and since <r(M, z) oc Go(z), the 
abundance of massive clusters is exponentially sensitive 
to the growth function. 

Additional details about the cluster abundance can be 
found in Sec. IVIDl where we describe the dependence of 
curvature estimates from SN, CMB, and cluster data on 
the modeling of the cluster growth information. 



D. Expansion rate 

The Friedmann equation gives the expansion rate as 
determined by the evolution of the density of various 
components: 

h(z) = H \n m (i + z) 3 + n r (i + zf 



PDEO) , n /, , x 2 ] 1/2 

+ S2k(1 + zj 



Pcr.O 



(11) 



where pde is the dark energy density and f2 m , f2 r , and 
£Ik are the present matter density, radiation density, and 
effective curvature density, respectively, in units of the 
critical density. Here we generalize the dark energy phe- 
nomenology to allow arbitrary evolution of the dark en- 
ergy density. For a general time-dependent dark energy 
equation of state w(z), the dark energy density evolves 
as 



Pde(z) = /O C r,0^DEexp 



3 / dz 

o 



,l + w{z') 



1 + z' 



(12) 



where J7de = 1 — — — is the present fraction 
of dark energy. 

The data we consider for curvature forecasts only have 
the ability to constrain the detailed dark energy evolution 
at low redshifts (z < z max where z max ~ 1.5), so we only 
allow complete freedom in pde(z) at late times. In mod- 
els with w(z) — 1, the dark energy fraction decreases 
rapidly with increasing rcdshift and therefore the exact 
form of the high-redshift dark energy evolution is unim- 
portant; for example, for flat ACDM, ^de^*) ~ 10~ 9 . 
However, models have been proposed in which the dark 
energy remains a significant fraction of the total density 
even at high z, such as scalar field models that track 
the density of the dominant matter and radiation com- 
ponents at early times [4l|, |42|, [43| . To account for such 
possibilities, wc parametrize early dark energy using an 
effective constant equation of state 

w(z > Z max ) = Woo, (13) 

as m IMHH09I . At z > Zmax, the dark energy density is 



PDe(^) = PDE(Zmax) 



1 + Z 

1 H~ -^ma: 



3(l+^oo) 



(14) 



Although this may not be a realistic model for early 
dark energy, it should be sufficient to absorb at least 
small effects of early dark energy on curvature estimates. 
Dark energy that behaves like a cosmological constant 
at z > Zmax would have Woo = — 1 and tracking mod- 
els have Woo = 0. Other values of Woo may provide an 
effective description of other types of models, e.g. track- 
ing models that transition to w < at z > z max can 
be approximated by a constant equation of state in the 
range — 1 < Woo < 0. We examine how well this approach 
works in the context of different early dark energy mod- 
els (as well as models with massive neutrinos, which have 
effects that are indistinguishable below the neutrino free 
streaming scale from early dark energy in the observables 
considered here) in Sec. IVIBI 

Since Woo only specifies the redshift evolution of the 
dark energy density at z > z max but not its normaliza- 
tion, the density at some redshift is required to com- 
pletely describe the early dark energy model. It is 
convenient to use the fraction of dark energy at z ma x, 
^DE(-Zmax), for this purpose. One way to estimate this 
quantity from SN data is to differentiate Df(z) at z max , 
since from Eqs. ^ and (fTTj) we can obtain 



PDE(Zmax) 
Pcr,0 



- E w 

i/DE 



3(l+w») 



(15) 



where E„ 



H(z max )/H = H^(dz/dD f )\ 2 
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E. Data for forecasts 

In this section, we summarize the main data assump- 
tions we use for the curvature forecasts presented in 
Sec. EU These data include future Type la supernova 
and CMB observations as probes of distances, and X-ray 
cluster abundances as a probe of the linear growth his- 
tory. The characteristics of the form er set of distance 
data match those assumed in IMHH09I . 

The supernova sample is taken to match the planned 
redshift distribution for the Super Nova/ Acceleration 
Probe (SNAP) Q, which covers redshifts 0.1 < z < 1.7. 
In addition, we assume a sample of 300 low-redshift SNe 
at 0.03 < z < 0.1. The intrinsic SN magnitude disper- 
sion is taken to be cr s tat = 0.15, and the systematic error 
is modeled as cr sys = 0.02[(1 + z)/2.7) in redshift bins of 
width Az = 0.1. Then the uncertainty in relative dis- 
tances from SNe in a Az = 0.1 bin with N SNe is 

a lnHoD = O^lnlOyOv-VLt+^ys- (16) 

The CMB distance priors we use are modeled on the 
specifications for the recently-launched Planck satellite 
H|. As in|MHH0j, we describe the CMB data with a 2D 
Fisher matrix F CMB including the distance to recombina- 
tion, D(z^), and the physical matter density, fl m h 2 . The 
elements of the covariance matrix C CMB = (^cmb-j-i 
are C™ B = (0.0018) 2 , C™ B = (0.0011) 2 , and C™ B = 
-(0.0014) 2 , where x = ln[D(z*)/Mpc] and y = n m h 2 . 

In addition to the distance constraints from SN and 
CMB data, in some cases we will consider the effect of 
addition al priors. For the MCMC analysis, we use the 
priors of lMHH09l which correspond to constraints from 
currently available data, including an 11% Hq prior from 
HST Key Project data [46[k a 3.7% BAO measurement 
of D(z = 0.35) from SDSSQ, and a 2.5% upper limit on 
the fraction of early dark energy at recombination from 
the WMAP temperature angular power spectrum [47| . 
For the growth reconstruction forecasts, we also consider 
future priors including a 1% measurement of Hq [48[ and 
a 1% upper limit on the dark energy fraction at recom- 
bination from CMB data, OrjE^*)- A stronger limit of 
~ 0.2% on the early dark energy fraction may be ob- 
tainable with future observations of CMB lensing in the 
context of specific parametrizations of early dark energy 
[S l50T ] , but here we adopt a more conservative prior to 
allow for the possibility that the limit may weaken upon 
including more general early dark energy behavior. 

Each of these constraints supplementing the SNAP and 
Planck data is implemented as a Gaussian prior with 
mean equal to the value in the fiducial cosmology. We will 
see in Sec. IVII that for a variety of cosmological models, 
particularly those in which the deviations in the high- 
redshift density evolution from the concordance model 
are mild, these additional priors are unnecessary for cur- 
vature estimates from future distance and growth data 
sets. 

For growth forecasts, we consider constraints on the 



amplitude and redshift evolution of galaxy cluster abun- 
dances using observations from the proposed Interna- 
tional X-ray Observatory (IXO) 0. The IXO is pro- 
jected to obtain 1 — 2% measurements of the growth func- 
tion Go(z) in Az = 0.1 bins over < z < 2, assuming 
that distances and the expansion rate are effectively fixed 
by other data sets such as SNe. When combined with 
the expected 1% measurement by Planck of the ampli- 
tude of scalar fluctuations from the CMB power spectra, 
the X-ray cluster data should also provide a ~ 1 — 2% 
measurement of G(z) at z < 2. 

We approximate the cluster growth information from 
IXO (and Planck for comparison of the cluster abun- 
dance and CMB power spectrum amplitudes) by a Gaus- 
sian likelihood with uniform uncertainties in Gq(z) at 
z < z raax in Az = 0.1 bins and G(z max ), using z max = 1.5 
corresponding to the maximum redshift at which dis- 
tances are well constrained by SNe. Our forecasts in 
Sec. IVII assume optimistic 1% growth uncertainties from 
IXO as the default assumption, but in Sec. I VI Dl we also 
examine how curvature constraints weaken for more pes- 
simistic assumptions and consider how degeneracies be- 
tween distances and growth in the cluster observables 
may affect the predicted curvature constraints. 



III. THE GEOMETRIC DEGENERACY 

In this section, we present examples of the degener- 
acy between curvature and dark energy evolution in the 
distancc-rcdshift relation, and show how growth infor- 
mation can break this degeneracy. 

The geometric degeneracy can take a variety of forms, 
depending on the dark energy modeling and the available 
data. For example, even for a cosmological constant A 
there is a degeneracy between Oa and Ok if the only 
input is the distance to recombination from the CMB 
(e.g., 0)- Adding more data breaks the degeneracy 
for the cosmological constant model, but the degeneracy 
persists for more complex dark energy models. Taken to 
an extreme, even if one has exact measurements of the 
distancc-rcdshift relation over the entire history of the 
universe, there is still a degeneracy if we allow arbitrary 
evolution of the dark energy density with redshift. This 
general form of the geometric degeneracy is what we focus 
on here. 

The degeneracy in distance data is apparent if we dif- 
ferentiate Eq. ([31) and use Eq. ((TTJ) to solve for the dark 
energy density: 

Pde(^) _ l + n K [g fid fl fid (*)]| N v 3(1+Wi) 

(17) 

assuming that SN data provide relative distance mea- 
sures [iJQ d D fid (z)]sN for some fiducial cosmology. [Note 
that this generalizes Eq. (fTS"]) .] Thus for any value of the 
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spatial curvature, there exists some dark energy evolu- 
tion that matches a given set of distance measurements. 

The ability of dark energy to match distances for any 
value of f2x is weakened if we introduce some mild restric- 
tions on the dark energy evolution, for example requiring 
that the dark energy density be nonnegative. Observa- 
tional constraints beyond the distance-redshift relation 
can also limit the possible dark energy behavior; for ex- 
ample, a dark energy fraction of a few percent or more 
at recombination would distort the CMB temperature 
power spectrum in ways that are not observed in WMAP 
data |,|47|,[I2. 

In practice, however, the redshift coverage of data is 
limited, so degenerate cosmological models only need to 
match distances at the redshifts where we can actually 
measure distances. Within the redshift "gaps" in our ob- 
servations, the dark energy evolution can deviate from 
Eq. Q17[) and still satisfy all available observational con- 
straints. 

Figure [5] shows the dark energy fraction, distance, and 
growth of four example models chosen to illustrate these 
points. These models all have nonzero curvature (two 
open and two closed), but fi m , h, and Woo are adjusted 
so that they all match relative SN distances HqD(z) at 
z < Zmax with z max = 1.7 (see middle panel of Fig. [2]) and 
CMB constraints on D(z*) and Sl m /i 2 for a flat ACDM 
model with fi m = 0.24 and h = 0.73. 

For each model, a dotted line shows how the dark en- 
ergy evolution would extend to z > z max if we required 
the models to match HqD(z) at all distances in the range 
< z < z* (but not the CMB absolute distance); for the 
closed models this would require a negative dark energy 
density, Q,-a^(z) < 0, at high z, and for the open mod- 
els J7de(^*) ~ 10 — 20% which violates current CMB 
constraints. However, due to the lack of observational 
constraints at z max < z < z*, not only can we find mod- 
els with nonzero curvature that match the SN and CMB 
data of the fiducial flat ACDM model, but there is in fact 
a set of degenerate models for each value of f^K with dif- 
ferent combinations of the matter density and early dark 
energy parameters. 

Specifically, to match an observed distance to recom- 
bination [D fid (z*)]cMB as well as SN distances, the con- 
straint that these models must satisfy [in addition to 
Eq. (T7|) at z < z max ] is 

- S^{KH^[Hi d D Rd (z^ SN }. (18) 

The bottom panel of Fig. [3] shows that despite being 
degenerate in the SN and CMB distance data, these mod- 
els are distinct in their growth evolution for different val- 
ues of the spatial curvature. This separation based on 
curvature is mostly independent of the {f2 m , u>oo} values. 
This indicates that growth observations can break the 
degeneracy between Q^, ^mi an d that is present in 
SN and CMB data and provide model-independent infor- 
mation about the curvature. 
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FIG. 2: Geometric distance degeneracy between curvature 
and dark energy dynamics. Curves show the dark energy 
fraction (top panel), relative distances {middle panel), and 
growth relative to z = (bottom panel) for models with 
distance-redshift relations that are degenerate with a fidu- 
cial flat ACDM model with f2 m = 0.24 and h = 0.73. Thick 
red curves correspond to open models with Qk = 0.01, 
and thin blue curves are closed models with Ok = —0.01. 
The solid and dashed curves only match the relative dis- 
tances at z < 1.7 (shaded region) and the CMB distance at 
z ~ 1000, but not distances at intermediate redshifts where 
the dark energy equation of state is assumed to be constant, 
w(1.7 < z < z*) = Woo. Dotted curves in the top panel show 
the dark energy evolution that would be required to match 
relative distances of the fiducial model at all redshifts. Pa- 
rameters of each model are given in the top panel legend; for 
each model, h is set so that Q m h 2 is the same as in the fiducial 
model. 



Note that a precise independent measurement of Hq 
(or, equivalcntly, combinations of SN relative distances 
with absolute distances, e.g. from BAO) determines the 
value of r2 m since the CMB precisely constrains f2 m /i 2 . 
This removes one parameter from the degeneracy de- 
scribed above, so the remaining degeneracy is between 
curvature and dark energy only with the matter density 
fixed. The freedom to match observations for any value 
of fix is thereby reduced, leading to constraints on cur- 
vature even in the absence of growth information. We 
discuss these constraints in relation to the curvature es- 
timates from distances and growth in Sec. IVIl 
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IV. MCMC METHOD 



For both the MCMC analysis presented in this sec- 
tion and the growth reconstruction method (Scc.fV]), it is 
convenient to view the constraints on curvature from dis- 
tances and growth in the following way. First, we assume 
that the distance-redshift relation is precisely measured 
at low z by SNe and at high z by the CMB. Given this 
distance data and assuming the validity of GR, we can 
use either the MCMC likelihood analysis or the analytic 
growth reconstruction technique to compute a predicted 
growth history that is consistent with the measured dis- 
tances. The main sources of uncertainty in the relation 
between measured distances and predicted growth are 
curvature and early dark energy (or massive neutrinos), 
with curvature primarily affecting the growth relative to 
z = [Go( z )] and early d ark energ y affecting the growth 
relative to high z [G{z)\ (|MHH09T ). Because of this de- 
pendence, measurements of Gq(z) yield constraints on 
curvature when they are compared with the predicted 
growth history, and measurements of G(z) help constrain 
deviations from matter domination at high z. 

Computing the distances and growth functions for 
each MCMC sample of the cosmological parameter space 
is straightforward as it only requires using Eqs. ([2|) 
and ([5]) to obtain the distance-redshift relation and solv- 
ing Eq. (|6|) for the growth. We will see in Sec.[V]that the 
growth reconstruction scheme is somewhat more compli- 
cated to implement. 

One of the main advantages of using the MCMC ap- 
proach is that the estimation of parameter uncertain- 
ties is also straightforward: as long as the MCMC sam- 
ples have converged to a stationary distribution approx- 
imating the joint posterior probability of the parame- 
ters (Sec. HVAj) . the marginalized probability for curva- 
ture can be obtained by binning the samples based on 
the value of fix- However, the parameter chains some- 
times converge quite slowly, especially for very general 
paramctrizations of dark energy combined with curva- 
ture. External priors like those described in Scc. lIIEl can 
improve MCMC convergence in some cases. 

A drawback of the MCMC analysis is that although we 
are trying to obtain model-independent results, we must 
still specify some model for the dark energy evolution. To 
provide a general parametrization at low z, we use sev- 
eral principal components of the dark energy equation 
of state as described in the following section. However, 
there is still some degree of unavoidable dependence on 
dark energy priors. In the absence of strong growth con- 
straints, the choice of priors can influence the curvature 
constraints as we will see in Sec. I VII The growth recon- 
struction method of Sec. [V] does not require specifying 
a parametrization for the dark energy equation of state 
and therefore suffers less from such effects. 



A. Growth predictions from distances 

The MCMC method of predicting the linear growth 
of perturbations is based on computing the growth func- 
tions of a large sample of cosmological models that fit 
observed distances well. The range of growth functions 
spanned by these models constitutes a prediction for the 
growth evolution base d on dista nces. This procedure is 
described in detail by IMHH09I . and here we present a 
summary. 

Two main ingredients are required: a parametrization 
of cosmological models and a description of how well 
these models fit the observed distances. Since we are 
concerned about possible degeneracies between dark en- 
ergy and curvature, we want to allow a wide variety of 
dark energy behavior. This is accomplished using a basis 
of principal components (PCs) for the dark energy equa- 
tion of state w(z) to parametrize general dark energy 
evolution below a redshift z m3 . x : 

w(z) - w Rd (z) = ^ diei(z), (19) 

i 

where the PCs are ordered according to how well they 
are measured by a particular combination of data. The 
PCs are eigenvectors of th e Fisher m atrix for the distance 
data, taken here and in IMHH09I to be the SNAP SN 
sample, CMB distance data from Planck, and current 
priors on Hq, D(z = 0.35), and Ode(z*) as described in 
Sec. Ill El Figure [3] shows the redshift dependence of the 
15 lowest-variance PCs. 




0.5 1 1.5 

z 

FIG. 3: The first 15 principal components of w(z) for SNAP 
and Planck (increasing variance from bottom to top), with 500 
redshift bins between z = and z m ax = 1-7 and uifid = — 1. 
The vertical dashed line shows the minimum redshift of the 
data assumed for computing the PCs, z^ n = 0.03. The PCs 
are offset vertically from each other for clarity with dotted 
lines showing the zero point for each component. 
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The Fisher matrix for the PCs is computed at some 
fiducial dark energy model specified by the equation of 
state WQd(z), usually taken to be a cosmological constant. 
The redshift range of the PCs is the same as the range 
of the SN data, with maximum redshift z max = 1.7. The 
PCs arc normalized as 



i=l 



iV s ,p C 
3=1 



i(zj)] 2 = N Z , PC , 



(20) 



where N Zt pc is the number of redshift bins, so that the 
shapes of the PCs are roughly independent of the chosen 
bin width. 

The highest- variance principal components have a neg- 
ligible effect on observable distances and growth due to 
their rapid oscillation in redshift. We therefo re trun- 
catc the sum in Eq. (19]) at 15 PCs, found by IMHH09I 
to be a sufficient number for a complete representation 
of the effects of dark energy variation at z < z max on 
the distance and growth observables. The dark energy 
description is completed by specifying the high-redshift 
evolution through the constant effective equation of state 
w(z > z max ) = Woo. Besides varying and the w(z) 
PC amplitudes in the MCMC analysis, we also include 
^m, Hq (parametrized through the combination £! m /i 2 ), 
and f^K as MCMC parameters, so there are 19 parame- 
ters in all. 

Using top-hat priors on the PC amplitudes, one can 
restrict the value of w(z) to a particular range, conser- 
vatively erring on the side of allowing too many models 
rather than too few. The priors corresponding to the 
< w < Wmax are a[~^ < on < a[ , where 



range u> n 



„(±) - 



N z ,pc 



77T ^ [( w mm + W max - 2w fid )e i (z j ) 

m x pc ■ 

3 = 1 

±(Wmax - W min )\ei(Zj)\], 



(21) 

assuming constant Wf\d(z) (|MHH09l ). 2 For example, re- 
quiring — 1 < w < 1 for quintessence models places strong 
limits on the allowed PC amplitudes. In the interest of 
keeping the dark energy evolution as general as possible, 
he re we wi ll only consider the weakest priors on w used 
in lMHH09l corresponding to — 5 < w < 3. 

By varying the cosmological parameters, computing 
the likelihood of each model for the assumed distance 
data sets, and using the Metropolis-Hastings criterion 
for deciding whether or not to accept a proposed step 
in the parameter space, the resulting set of parameter 
combinations trace out the joint posterior probability of 
the parameters (e.g., [Hi], [H, HI]). To determine when 
the number of MCMC samples is large enough that the 
parameter chains have converged to the posterior distri- 
bution, we run 4 independent chains and require that the 
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2 Similar top-hat priors are derived in Ref. [53j in the context of 
principal components of the reionization history. 



01234 01234 
z z 



FIG. 4: Growth functions of MCMC samples with general 
dark energy equation of state variations at z < 1.7 within the 
range —5 < w < 3, including either early dark energy (EDE) 
at z > 1.7 (woo ^ —1; top), curvature (J1k 0; middle), 
or both (bottom). The left panels show growth relative to 
early times, and the right panels show growth relative to the 
present. Dashed red curves show growth in the fiducial flat 
ACDM model (fi m = 0.24, h = 0.73). Samples are selected 
randomly from those with likelihoods satisfying A\ 2 < 4, 
but for visual clarity we plot samples that are approximately 
evenly spaced in G(z — 0) (left) or Go(z = 4) (right). The 
dotted vertical line in each panel marks the division between 
the low-z and high-z dark energy descriptions at z = 1.7. 



Gelman-Rubin statistic satisfy R — 1 < 0.01, indicating 
that the variance of the mean value of a parameter be- 
tween different chains is much smaller than the variance 
within a single chain (57| . 

This MCMC procedure produces a variety of cosmolog- 
ical models that fit the fiducial distance data reasonably 
well. For each of these models, we compute the growth 
history using Eq. ([6]). The distribution of the resulting 
growth functions then forms the prediction for growth 
from distances. 

Plotting sets of predicted growth functions from 
MCMC analyses with different degrees of freedom in 
the cosmological models reveals how distance-matched 
growth functions depend on curvature and early dark 
energy. Figure 2] shows the growth evolution of selected 
MCMC samples in chains with early dark energy vary- 
ing and curvature fixed to fix = 0, curvature varying 
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and early dark energy fixed to Woo = — 1, or variation 
in both curvature and early dark energy (in addition to 
variation in the low-z dark energy equation of state via 
principal components). Early dark energy mainly affects 
the growth amplitude relative to high redshift, G(z), with 
very little effect on the shape of the growth evolution at 
low z characterized by G$(z). On the other hand, curva- 
ture strongly influences Gq(z) but has less of an effect on 
G(z): for a fixed distance-redshift relation, open (closed) 
models have larger (smaller) Gq(z) than a flat universe. 
This is the same effect that we see in the geometric de- 
generacy examples in Fig. [21 

This difference between the effects of curvature and 
early dark energy supports the claim that measurements 
of Go (z), combined with distance measurements, can 
constrain curvature with little dependence on early dark 
energy or other high-rcdshift phenomena. Measurements 
of G(z) can be used to place limits on any residual effects 
that early dark energy might have on these curvature es- 
timates. 

One interesting feature of the spread of the distance- 
matched Go(z) evolution due to curvature is that it tells 
us how precise growth measurements must be to improve 
on model-independent curvature constraints. In particu- 
lar, for the forecasts shown here, the uncertainty in the 
predicted growth at z ~ 1 (at 68% CL) is about 5 — 10% 
(|MHH09l ). 3 so any direct measurements of the growth 
evolution that are less precise than this will not appre- 
ciably reduce the uncertainty in J7k- For example, cur- 
rent cluster measurements of growth are not yet precise 
enough to significantly reduce the uncertainty in curva- 
ture from current distance constraints (40| . 

As observed in IMHH09I . MCMC estimates of fi K from 
distance data are dependent on w{z) priors if "phantom" 
dark energy models (w < — 1) are allowed. This depen- 
dence is essentially a volume effect related to the large 
volume of degenerate w < — 1 models that are correlated 
with Qk < [58| . Precise growth measurements help 
reduce this dependence on dark energy priors as we will 
see in Sec. IvTl 



B. MCMC estimates of curvature 

The procedure described in the previous section pro- 
duces chains of parameter combinations that match the 
SN and CMB distance data, as well as the additional pri- 
ors. To obtain forecasts for J7k from distance and growth 
data, we need to add the growth information from mea- 
sured cluster abundances. 

The first step for including the growth constraints is 



to simply compute the growth evolution for each MCMC 
sample. With a growth history associated with each 
distance-matched MCMC sample, we then use impor- 
tance sampling of the parameter chains [59| to reweight 
the samples according to the growth likelihood described 
in Sec. Ill El The posterior probability for fix from the 
distance and growth forecasts is then computed as usual 
by marginalizing over the other MCMC parameters. In 
Sec. I VII we will describe the resulting curvature con- 
straints and compare them with the forecasts from the 
growth reconstruction method. 



V. GROWTH RECONSTRUCTION 

The MCMC approach in the previous section used a 
general parametrization of the dark energy equation of 
state to predict growth evolution from measured dis- 
tances by searching for w(z) satisfying the distance con- 
straints and computing the corresponding growth evo- 
lution. However, it is possible to go from distances to 
growth directly without using an intermediary like w(z). 
Here we will first summarize this growth reconstruction 
method and then show how it can be used to study 
model-independent curvature constraints. Additional de- 
tails about the procedure we use for growth reconstruc- 
tion from simulated SN and CMB data are provided in 
Appendix [B] 

The growth reconstruction equations derived by Alam, 
Sahni, and Starobinsky [23l.[60| express the amplitude of 
linear perturbations 5 as a function of the comoving dis- 
tance assuming spatial flatness, Df [Eq. (J5J)]. For con- 
venience, wc define a dimcnsionlcss comoving distance, 
X = H Df. 4 

Starting from Eq. ([5]) for the linear growth function, 
we can rewrite the equation in terms of x by using d/dt — 
—Hq(1 + z)d/d\ and writing the matter density as p m = 
3r2 mJ ff J /(87rG A r)(l + zf to get 

H^(l+z)[(l+z)8']'-2H H(l+z)S / -^H^n m (l+z) 3 S = 0, 

(22) 

where primes denote derivatives with respect to x- Since 
dx/dz = Hq/H, we can replace H by H dz/dx = H (l + 
z)' . After dividing by Hq(1 + z) 3 , this yields 

{(I + z)5>}>- 2(1 + zJl = 3 n s 



(l + zf 



(23) 



The left hand side of this equation is equal to [(1 + z)8 h 
(1 + z)'5']/(l + z) 2 , which can be written as 



(24) 



3 Note that the uncertainty in growth predictions depends on the 
priors assumed in addition to the SN and CMB data forecasts, 
in this case taken from current Ho, BAO, and CMB data as 
described in Sec. Ill El 



4 Note that in Ref. \2.j this quantity is called E, but here we 
use the notation \ instead to avoid confusion with the common 
definition E(z) = H(z)/H . 
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Integrating Eq. (|24[) over \ leads to the main growth re- 
construction equation p3l . |60| , which is an integral equa- 
tion for <5(x): 



<Kx) = l + «o / dXi[l + «(Xi)] (25) 
Jo 

/ dXl[l + ^(Xl)] / d*2 «(X2), 

/o Jo 



+ 2^ 



where 6' is the derivative of 5(x) at z = 0, and the growth 
function is normalized to <5(z = 0) = 1. Solving for the 
growth function involves making an initial guess for 5(x) 
and plugging it into the right hand side of Eq. ([23]) , taking 
the resulting 6(x) from the left hand side and plugging 
it back in to the right hand side, and repeating until the 
solution has converged. 

As noted in Ref. [23[, Eq. (|25|) only requires integra- 
tion of observed quantities, so the growth reconstruction 
method is more stable to the presence of scatter in the 
data than other methods that use derivatives of the ob- 
served distance-redshift relation. Differentiation uses in- 
formation from the data over only a small range of red- 
shift or Xi whereas an integral from to x uses data over 
the entire range which helps average out statistical noise. 
For x near the number of data points used is still small, 
but the reduced statistical power can be offset by having 
smaller intrinsic uncertainties at lower rcdshifts. 

By differentiating Eq. (j2l)j) and evaluating it at x = 
x(zmax), we can shift the boundary condition on the 
derivative of S(x) from z = to z = z max : 



6' 

c7 max 

Q= l + z~ 



where <5^ ax = 5'(x(z E 



2 



X(*max) 



d X S(x), (26) 



ten as 



<)). Then Eq. ([25]) can be rewrit- 



1 



d X i[l + z{xi)} 



3 rx rx( z n 
-rj m y dxi[l + z(xi)] 



(27) 
dX2 S(X2)- 



This is the form of the growth reconstruction equation 
that we will use for the curvature forecasts. 

One advantage of setting the boundary condition for 
<5'(x) at x( z max) instead of x = is that S' max de- 
pends mainly on the assumed cosmology at high red- 
shifts, whereas 5' depends on not only the high-z as- 
sumptions but also the low-z SN data constraints and 
iterative growth solution. Also, setting 5' max makes it 
easier to ensure that S(x) is smooth at x(z max ) by re- 
quiring the same derivative there for both the fiducial 
model at z > z max and the reconstructed growth func- 
tion at z < z max . Rather than setting the value of the 
5'(x) boundary condition using an integral over S(x) as 
suggested in Ref. [60l |. here we use an approximate an- 
alytic form for the high-redshift growth function valid 
when matter is the dominant component. 



Solving Eq. (|27|) for 8(x) requires specifying the func- 
tion z(x) and three parameters: fix, ^m, and <5^ ax . Some 
of these inputs to the growth reconstruction are set by 
the SN and CMB data constraints, while others are free 
parameters. Because of this link between the distance 
data sets and the growth reconstruction, some of the scat- 
ter in the SN and CMB observations will propagate to 
uncertainties in the reconstructed growth. We compute 
the mean values of reconstructed growth observables and 
their uncertainties using Monte Carlo simulations of the 
SN and CMB data, modeled on SNAP and Planck as for 
the MCMC method. 

The following steps summarize the procedure for these 
Monte Carlo simulations. Additional details about each 
step are provided in Appendix [Bl 

1. Assume values for the curvature and any additional 
parameters describing the high-redshift cosmology; 
here we use Woo for parametrizing early dark en- 
ergy. The output of the Monte Carlo simulations 
will be the conditional probability for the growth 
observables given fix and u>oo- 

2. Draw a realization of the SN data \HqD(z)\ for 
the fiducial cosmology and estimate z(x) from the 
data. Given an assumed value of Qk, we can invert 
Eq. © to get 



1 



.Q-l 



(28) 



where (x) is sinh 1 x for an open universe, 
sin -1 x for a closed universe, and x if Oj< = 0. To 
reduce bias in the estimated x( z ) relation, we take 
the maximum redshift for the growth reconstruc- 
tion to be z max = 1.5, slightly lower than the max- 
imum SN redshift of z = 1.7 that was used as z max 
in the MCMC method of Sec. El (see Appendix [0 
for details). 

3. Draw a realization of the CMB data [D(z*) and 
J! m /i 2 ] and use Eq. (fTg|) along with the assumed 
values of fix and Woo to compute the value of Cl m 
required by the CMB parameters. 

4. Steps 1 — 3 fix the cosmological model at z > z max , 
so we can compute G{z max ) for this cosmology as 
one of the growth observables, and use the ap- 
proximate high-z growth solution of Appendix [A] 
[Eq. (|A5|) ] to set S' max for the growth reconstruc- 
tion. 



5. Solve Eq. (|27|) to find the reconstructed growth, 
S(x), for the particular realization of SN and CMB 
data. Use the x( z ) relation from step 2 to express 
this solution as 5(z) or Gq(z). 

Repeating these steps for many realizations of the dis- 
tance data produces a distribution of the growth observ- 
ables g, which include G(z max ) from step 4 and Gq(z) 
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FIG. 5: Effect of early dark energy on distance and growth. 
Top panel: Difference between the distance to recombination 
and the distance to z mBX = 1.5. Bottom panel: Growth at 
Zmax relative to recombination. Early dark energy models 
are parametrized by the fraction of dark energy at z ma x and 
Woo = w(z > z max ). From top to bottom in each panel, 
ujoc = -1, -0.7, -0.5, -0.3, -0.2, -0.1, and 0. Flat ACDM 
is assumed here with Q m — 0.24 and h = 0.73. 

from step 5, for the chosen values of £Ik and w^. This 
procedure can be carried out at several different values 
of £Ik and Woo to map out the conditional probability 
P(g|£lK, Woo). This probability describes the growth pre- 
dictions from distance data in the context of the growth 
reconstruction method. 

Note that in step 3, there will generally be degenera- 
cies between fix, S^m, and early dark energy in the CMB 
constraints, increasing the uncertainty in the modcl- 
indcpcndent estimate of curvature. Fortunately, mea- 
surements of growth relative to high z can reduce this 
uncertainty. The top panel of Figure [5] shows the effect 
of various early dark energy models on the difference be- 
tween the distance to recombination and the distance to 
Zmax- Both Woo > -1 and r2 DE (z max ) > 0.05 are re- 
quired for the distances to be significantly affected. As 
the bottom panel of Fig. O shows, the growth relative to 
high redshift, G(zmax), is sensitive to early dark energy 
in a way that is similar to how the distances depend on 
early dark energy. Therefore, precise measurement of the 
growth at z max relative to growth at early times (for ex- 
ample, by comparing the normalization of growth from 
cluster abundances with the amplitude of CMB power 
spectra) can constrain the effect of early dark energy on 
high-z distances, thereby reducing the uncertainty in the 
^m _ ^K relation from the CMB constraints. Although 
G(z) predicted from distances is relatively insensitive to 




FIG. 6: Contours of the growth rate / at z m ax = 1.5 as a 
function of fix and fiDB(zmax), for three choices of uioo at fixed 
fi m =0.24 {top) and three choices of fi m at fixed w x = — 1 
[bottom). 



curvature in a direct sense (see Fig. [5]), it provides an im- 
portant complementary constraint to the curvature esti- 
mates from distances and Gq(z) that reduces model de- 
pendence. 

The approximation for <5 max used in step 4, which 
is described in Appendix [Aj assumes that the curva- 
ture and early dark energy fractions at high z are small 
enough that the growth can be written as a perturba- 
tion to the matter-dominated solution where G(z) is con- 
stant. The dependence of 6' max on early dark energy 
turns out to be fairly weak. Figure [5] shows the value 
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FIG. 7: Growth reconstruction for S1k = and Woa = — 1 
where the fiducial cosmology is flat ACDM with Q m = 0.24 
and h = 0.73. The shaded band is the 68% CL region for 
Go (2) reconstructed from simulated SNAP and Planck data, 
and the true growth evolution is plotted as a solid curve. The 
redshift bin width for the reconstruction is Az = 0.05. 



of the growth rate / = d In S/d In a = 1 + d In G/d In a at 
^max aS a function of curvature and the fraction of dark 
energy at z max . Changing Woo only slightly shifts the 
contours of /, especially for models that satisfy CMB 
constraints on the dark energy fraction at recombina- 
tion. Therefore, at fixed curvature S' max /5(z Taax ) = —(1 + 
Zmax) _1 £(zm a x)/(zma) mainly depends on S! D E(Zmax), 
which can be estimated from the SN data using Eq. ([15)) . 
Likewise, for curvature at the level of a few percent or 
less, uncertainty in f2 m does not strongly affect ^ax- 
In summary, the growth reconstruction method as out- 
lined above requires choosing values for two parameters: 
the curvature S1k and the early dark energy equation of 
state Wqo. All other quantities are estimated from the dis- 
tance information provided by measurements of SNe and 
the CMB. Given this distance data, the end result of the 
growth reconstruction procedure is a prediction for the 
growth evolution for each {fix, Woe} pair, P^g^x, w^). 

Figure [7] shows an example of the reconstructed growth 
evolution from the distance data for the fiducial flat 
ACDM model, assuming fix = and w^, = — 1. In 
the lowest redshift bins, the growth reconstruction is bi- 
ased due to the relatively small number of SNe in the 
SNAP distribution at low z, but over most of the red- 
shift range the reconstruction is in good agreement with 
the true growth function. 



A. Growth reconstruction estimates of curvature 

By combining the predicted growth from distances 
Pi (g|^K) with measurements of the growth observables 
g, we can obtain an estimate of Qk from the growth re- 
construction method: 5 



p(o K )= / dgP d (n K , g )p g ( s ), 



(29) 



where P d (fi K ,g) = Pz(g|^K)P P rior(^K) is the joint prob- 
ability for curvature and growth observables from the 
growth reconstruction from distance data as described in 
the previous section combined with any additional prior 
information about SIk, and P g (g) represents the con- 
straints on the growth observables from probes of the 
growth history such as clusters. We use subscripts d 
and g to indicate constraints coming from distance and 
growth data, respectively. 

Assuming that the Monte Carlo simulations produce 
conditional probabilities for growth observables at fixed 
Ok that can be approximated as a multivariate Gaussian, 
we can write 



Pd{g\n K ) 



1 



(2 7 r)"/ 2 (detF d )- 1 /2 



exp ( -i<5g T F d (5g 



(30) 

where Fd is the Fisher matrix for the n growth observ- 
ables from the simulated distance data, computed by in- 
verting the covariance matrix from the growth recon- 
struction Monte Carlo simulations, and <5g = g — gd is 
the deviation of the growth observables from the mean 
growth reconstruction solution. Note that both F^ and 
grf depend on the value of I^k (and additional parameters 
such as Woo). 

If the growth observations are also well approximated 
by a multivariate Gaussian with covariance matrix C ff , 



Fisher matrix F, 



r 1 



and mean values g 9 (assumed 



to be equal to the true growth history of the fiducial 
model), then the posterior probability for fix in Eq. ([29]) 

is 



P(fi K ) oc P pri or(ft K ) [detCE + F^Fg)]- 1 / 2 



(31) 



x exp 



i A T F 9 (Fd + F g )~ l F g A - ± A T F 9 A 



where A = g^ — g ff is the difference between the average 
growth evolution predicted from distances at a particular 
assumed value of £Ik and the true growth evolution, and 
I is the n x n identity matrix. We describe the resulting 
forecasts for curvature in the next section. 



5 For notational compactness, in this section we suppress depen- 
dence on additional parameters varied in the growth reconstruc- 
tion analysis such as u>oo, but in general, these should appear 
along with S1k wherever there is dependence on curvature. 
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VI. MODEL-INDEPENDENT CURVATURE 
CONSTRAINTS 

In this section, we use the techniques for combin- 
ing distance and growth measurements described in Sec- 
tions QV] and [V] to obtain forecasts for spatial curvature 
constraints from the simulated SN, CMB, and cluster 
data of Sec. UTEl 

The accuracy of curvature estimates depends not only 
on the assumed characteristics of the distance and growth 
data sets, but also on the fiducial, "true" cosmological 
model assumed for the forecasts. We begin in Sec. IVI Al 
with the simple case in which the fiducial model is flat 
ACDM, and use this example to compare forecasts from 
the MCMC and growth reconstruction methods. In 
Sec. IVI Bl we generalize to other fiducial cosmologies that 
are more or less consistent with constraints from current 
data. Tests of the dependence of these results on model- 
ing of the SN, CMB, and cluster data sets are described 
in Sections IVI CI and IVIDI Finally, Sec. IVI El compares 
these distance and growth constraints on curvature with 
other model-independent tests of curvature. 



A. Flat ACDM 

We begin by comparing MCMC and growth recon- 
struction constraints on f?K hi the context of a flat 
ACDM fiducial cosmology, taking the parameters to be 
51 m = 0.24 and h = 0.73. Figure [8J shows curvature fore- 
casts for two different sets of data. Both cases include 
distance constraints from SNe modeled after SNAP and 
CMB data based on Planck. On their own, these data 
sets would place almost no limits on curvature, assum- 
ing that general forms of the dark energy evolution are 
allowed. 

In the top panel of Fig. [SJ the SN and CMB data are 
supplemented by weak priors based on current data: an 
11% H prior from HST Key Project data [Ilk a 3.7% 
BAO measurement of D(z = 0.35) from SDSS 2|], and a 
2.5% upper limit on the fraction of early dark energy at 
recombination from the WMAP temperature spectrum 
[47l ]. We adopt these priors here for consistency with 
the MCM C analysi s of model-independent growth pre- 
dictions in IMHH09I . 

Using only the current priors in addition to the SN and 
CMB data, curvature is determined with an accuracy of 
a(flK) ~ 0.02. However, the results depend on the analy- 
sis method. In particular, the MCMC constraints on 17k 
depend stro ngly on th e assumed priors on the dark energy 
parameters (|MHH09f ). For priors that are flat in the w{z) 
principal component amplitudes [P(oti) = constant], the 
curvature constraint is biased toward J7k < 0. Taking 
alternate priors which are instead flat in the density as- 
sociated with each principal component at z max = 1.7, 
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FIG. 8: Comparison of MCMC and growth reconstruction 
curvature forecasts for the flat ACDM fiducial cosmology. 
Top panel: combining current priors on Ho, D(z = 0.35), 
and Qt>e(z*) with SN and CMB forecasts based on SNAP 
and Planck, respectively. Bottom panel: combining SNAP 
and Planck forecasts with cluster forecasts based on IXO as 
a probe of the growth evolution. For comparison, the dotted 
curve shows a forecast for SNAP, Planck, and future 1% mea- 
surements of Ho and Que(z,) (without growth constraints). 



the resulting curvature constraint is less biased, but has 
a long tail at f2x > 0. The curvature estimate from 
the growth reconstruction method using the same data 
is more consistent with the flat-density prior constraint 
from the MCMC analysis, but has a slightly narrower 
distribution. 

The lower panel of Fig. [5] shows the resulting curva- 
ture forecasts when we drop the current priors on Hq, 
BAO distance, and the early dark energy fraction, and 
instead combine the SN and CMB distances with 1% 
growth function measurements from IXO clusters. Us- 
ing growth to break the geometric degeneracy in dis- 
tances, the curvature constraint improves to an accuracy 
of (j(£!k) = 0.0022. Not only is this a much stronger 
constraint on curvature than without the growth infor- 
mation, but it is also significantly less dependent on the 
parameter estimation methodology. The MCMC con- 
straints using different types of dark energy priors are 
consistent with each other and unbiased, and the fore- 
cast from growth reconstruction agrees with the MCMC 
results. 

Since we are taking a forecast for future data as the 
source of growth information, it is not quite fair to com- 
pare this future distance plus growth constraint with the 
curvature estimate from current Hubble constant, BAO, 
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FIG. 9: Forecasts in the Qk-w^, plane for the flat ACDM 
fiducial model, including various combinations of growth in- 
formation in addition to SN and CMB distance constraints: 
G(zmax) (dashed contours), Gq(z) at z < 1.5 (solid con- 
tours), or both (shaded blue contours). The combined 
growth constraint without the covariance between the pre- 
dicted G(z m ax) and Go(z) from the growth reconstruction is 
shown as shaded red contours. Contours are plotted at 68% 
CL (thick curves/dark shading) and 95% CL (thin curves/light 
shading). 



and early dark energy priors. We use those priors only 
to illustrate the possible dependence of estimates that 
do not use growth information on dark energy priors or 
other details of the analysis, and for comparison with the 
earlier MCMC results from IMHH09I . If we instead com- 
bine SNAP and Planck forecasts with future priors on 
H and Ode (2*)) both with 1% accuracy (see Sec. Ill E|) . 
the curvature estimate improves but remains weaker and 
more skewed than the distance plus growth forecasts (see 
bottom panel of Fig. [5]). We discuss the impact of future 
BAO measurements and other types of data on modcl- 
independent curvature constraints in comparison to con- 
straints from distances and growth in Sec. IVI El 

The MCMC and growth reconstruction forecasts show 
that combinations of distances and growth have the po- 
tential to provide reliable, model-independent measure- 
ments of curvature with ~ 0.2% accuracy, at least in the 
context of a flat ACDM cosmology. Before exploring the 
forecasts for this method for other fiducial cosmological 
models, we examine how different parts of the distance 
and growth data contribute to the curvature estimate. 

Figure [5] shows the impact of different components of 
the assumed growth information from IXO cluster data 
on the joint constraint on curvature and early dark en- 



ergy. Recall from Sec. Ill El that there are two types of 
growth information used here: 1% measurements of the 
growth evolution at < z < z max that probe Gq(z), and 
a 1% measurement of the growth at z max relative to the 
growth at recombination, G(z max ). The unshaded con- 
tours in Fig. [5] show the results of using each of these 
growth constraints separately (in addition to the SN and 
CMB data). 

If the growth information consists only of Gq(z) (i.e. 
just the relative evolution of the cluster mass function at 
low z with unknown normalization relative to the CMB), 
the constraints on f2x are a few times weaker than the 
combined growth constraints for < — 1. As Wrx, ap- 
proaches 0, the Go constraint shifts toward open models. 
As a result, P(f2jc) marginalized over Woo is weakened 
even further and is skewed toward J7k > 0. 

On the other hand, growth information from G(z max ) 
only (e.g., comparing o~s from clusters and the CMB) 
places very weak constraints on curvature, only signifi- 
cantly limiting the range of allowed open models. How- 
ever, the information provided by G(z max ) is complemen- 
tary to that from Go(z) since it cuts off the Go degen- 
eracy at fix > and Woo > —0.5. By removing this 
degeneracy, the combined growth constraints are much 
less sensitive to the early dark energy parameters, result- 
ing in a stronger, unbiased curvature estimate. 

Moreover, the covariance between G(z max ) and Go (z) 
predictions from measured distances reduces the curva- 
ture uncertainty beyond what would be expected from 
simply combining the separate G(z max ) and Gq(z) con- 
straints. As shown in Fig. [51 the distance plus growth 
constraints on curvature with the G(z max ) — Gq(z) co- 
variance removed by hand (red shaded contours) are a 
few times weaker than the full constraint including this 
covariance (blue shaded contours). 

The covariance between the growth observables comes 
primarily from the CMB distance constraint. Matching 
the distance to recombination requires a balance between 
the low-z and high-z dark energy evolution. For exam- 
ple, increasing the dark energy density at low redshifts 
tends to decrease D(z Sf ), and decreasing the high-redshift 
dark energy density can compensate for this shift. The 
CMB distance priors therefore anticorrelate the dark en- 
ergy evolution at early and late times. This results in a 
positive correlation between the predicted values of Gq{z) 
and G(z max ) from SN and CMB data. In the example 
above, the larger dark energy density at low z enhances 
the suppression of the late-time growth of perturbations, 
which corresponds to higher values of Gq(z). A smaller 
dark energy density at high z results in less growth sup- 
pression at early times and therefore a higher value of 
G(z max ) . 

This positive correlation between Go and G for growth 
reconstructed from distances is the opposite of the effect 
of curvature. Relative to a flat universe, an open ge- 
ometry tends to suppress growth at both early and late 
times while a closed geometry has the opposite effect. 
This means that increasing J7k increases Gq(z) while dc- 
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creasing G(z max ), and vice versa. The positive covariancc 
between these two types of growth observables required 
by CMB constraints therefore leads to stronger limits on 
curvature. 

To summarize, in the case of flat ACDM the primary 
information about curvature from growth combined with 
distance data comes from Go (z) , as expected based on the 
discussions in the previous sections. However, normaliza- 
tion of low-redshift growth relative to the high rcdshift 
of recombination can significantly improve the accuracy 
of the curvature estimate by reducing early dark energy 
uncertainties and through the covariance with the low-z 
growth evolution required to match SN and CMB dis- 
tances. 



B. Dependence on cosmology 

We now turn to curvature forecasts from distances 
and growth for fiducial cosmologies other than the flat 
ACDM example of the previous section. Since the an- 
alytic growth reconstruction method is more efficient at 
exploring different assumptions about the true cosmol- 
ogy and properties of the data, we will primarily rely on 
that method for the forecasts in this section. 

As the first test of dependence on the fiducial cosmol- 
ogy underlying the data, we consider flat ACDM mod- 
els with different parameters from those in the previous 
section. For a model with f2 m = 0.3 and h = 0.7, the 
constraint on curvature is almost identical to that for 
fl m = 0.24 and h = 0.73. Thus variation of flat ACDM 
parameters in the range allowed by current data does not 
significantly affect the curvature forecasts. 

Next, we consider models with different values of the 
curvature. We saw in the previous section that if the 
universe is actually flat, future distance and growth 
probes can exclude at 95% CL alternate models with 
| Ok | > 0.005. It is also interesting to ask whether true 
nonzero curvature could be detected (i.e. distinguished 
from fix = 0) using this method. Figure [TU] compares 
the forecast for the flat ACDM model of the previous 
section with forecasts for an open model and a closed 
model, both with \Qk\ = 0.01 and the other parameters 
unchanged. The resulting constraints on curvature for 
the open and closed models are nearly identical to that 
for the flat model, except for being shifted to be centered 
on the true curvature. Therefore, distances and growth 
from SNAP, Planck, and IXO data enable percent-level 
curvature to be cleanly distinguished from flatness. 

Curvature forecasts are similarly precise and unbiased 
for models with different dark energy evolution from 
ACDM at low redshifts. For example, if we adopt the 
commonly-used parametrization of the dark energy equa- 
tion of state w(z) = wq + w a z/(l + z) [6l|, [&2j with 
wo = —0.8 and w a = —0.5, and keep all other param- 
eters the same as in the fiducial fiat ACDM model, then 
the growth reconstruction analysis again returns an un- 
biased estimate of Qk with (j(£Ik) ~ 0.002. Note that 
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FIG. 10: Growth reconstruction forecast of the marginalized 
probability of Ok for models with Qk = (solid), Ok = —0.01 
(short dashed), and Ok = 0.01 (long dashed). Each model is 
a ACDM cosmology with n m = 0.24 and h = 0.73. 



in this model, dark energy is negligible at high redshifts 
since lim z _ >00 w(z) = wo + w a which is less than —1. 

Although changing parameter values in the context of 
flat ACDM and similar cosmological models has little ef- 
fect on the accuracy of the estimated curvature from fu- 
ture distances and growth, changing the fiducial model 
for the high-rcdshift universe can have more interesting 
consequences. We consider first early dark energy sce- 
narios and then models with massive neutrinos. In both 
cases, the G(z max ) constraint takes on a much larger role 
and constraints from Gq(z) alone become unreliable. 

Given the fact that current data are consistent with 
flat ACDM, models with significant early dark energy 
are more likely to be open than flat or closed since the 
geometric effect of negative spatial curvature can com- 
pensate for the reduced distance to recombination due 
to early dark energy (e.g., [1]). Figure QT] shows the 
SNAP+Planck+IXO constraints on a wq— w a model with 
wq + w a ~ to act as early dark energy and Qk = 0.025 
compared with the flat ACDM model of the previous sec- 
tion. 

For the open early dark energy model, the best fit to 
distance and growth constraints is near the true param- 
eter values (the effective value of w at z > z max that 
matches the fiducial CMB distance and growth function 
normalization is ss —0.25). However, a significant 
degeneracy between curvature and early dark energy re- 
mains even with the combination of distance and growth 
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FIG. 11: Forecast in the Qk-w x plane for an open universe 
with early dark energy: Qk = 0.025, Wo = —1.08, w a = 1.02 
{shaded blue contours, right). The mean value of G(z ma x) for 
the growth function reconstructed from the fiducial distances 
of this model is plotted with dashed, unshaded contours at 
G(zmax) = 0.5, 0.6, 0.7, 0.8, and 0.9, from top to bottom. The 
flat ACDM forecast from Fig. [11] is plotted for comparison 
[shaded gray contours, left). Crosses mark the best fit point 
for each model. Shaded contours are plotted at 68% and 95% 
CL. 



data. One consequence of this is that the curvature con- 
straints marginalized over w^, would be misleading; for 
example, the shape of the wq — w a model constraints 
in the £Ik — Woo plane cause P(Q,k) to be biased high, 
although maximizing the likelihood over Woo instead of 
marginalizing would reduce the apparent bias somewhat. 
Furthermore, the marginalized Qk constraints are weaker 
than the 0.2% fiat ACDM estimate, despite the fact that 
the contours in the £1k — Woo plane arc of comparable 
width. 

Despite the loss of precision in the curvature estimates 
in the presence of early dark energy, the constraints from 
distances and growth for such models remain extremely 
interesting since they provide clear evidence for either 
large amounts of early dark energy or nonzero curvature, 
if not both. And, in fact, for models like this example 
with J1k > 0.02 it is still possible to confidently exclude a 
flat universe independent of the assumptions about early 
dark energy. 

Unlike the models without early dark energy where 
the main curvature constraints come from distances plus 
Gq(z), and G(z max ) has a lesser role (see Fig. [5]), for 
early dark energy models the G(z max ) constraint com- 
bined with distance data is the main source of both the 



early dark energy and curvature constraints. In fact, 
curvature estimates using distances and measurements 
of Gq(z) alone would indicate a preference for flatness 
and no early dark energy even if the true model was the 
wo — w a example with f^K = 0.025. Other types of cur- 
vature constraints, for example using a 1% Hq prior in 
addition to SNAP and Planck data, similarly fail to de- 
tect the nonzero curvature and early dark energy for such 
models. While a measurement of G(z max ) is a helpful ad- 
ditional constraint on models in which dark energy only 
becomes important at late times, it is crucial for obtain- 
ing accurate constraints on models with early dark en- 
ergy. Figure [11] shows contours of G(z max ) for the open 
wq — w a model; the forecast roughly follows these con- 
tours, but is tilted somewhat due to the Gq(z) — G(z max ) 
covariance from the CMB distance constraint. 

To test how well the early dark energy parametrization 
with Woo can model other forms of early dark energy evo- 
lution, we use the parametrization of Ref. [63j | in which 
the dark energy fraction approaches a constant value O e 
at high rcdshift: 
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where wq = w(z = 0). (This form assumes f2j< = 0, so 
we only use it to simulate distance and growth data for 
a flat universe.) The equation of state at high redshifts 
during matter domination is w ~ 0. However, since the 
transition between w ~ wo at low z and w ~ at high z 
can occur at z ^ z max the effective value of Woo that best 
matches distance and growth observables for this model 
is not necessarily Woo = 0. 

We take the parameters for the simulated SNAP, 
Planck, and IXO data to be w = -0.9, fl e = 0.02, 
Sl m = 0.225, h = 0.73, and Q K = 0. As for the pre- 
vious early dark energy example, the constraint on Ok 
marginalized over Woo is biased due to the shape of the 
likelihood in the f2j< - Woo plane. However, for this model 
there is some additional bias such that the best fit is at 
JIk ~ 0.004. Given the increased width of the likelihood 
in f^K [cr(fiK) ~ 0.003] and the long tail toward 17k < 0, 
the true curvature f2x = is not strongly disfavored de- 
spite this bias. The bias in f2j< is likely a reflection of 
the limitations of the Woo parametrization of early dark 
energy, suggesting that a more flexible parametrization 
may be needed to accurately model a wide variety of early 
dark energy behavior. 

Additional priors can help reduce bias in curvature 
even with the default Woo modeling of early dark energy 
in the growth reconstruction analysis. The open models 
favored in the analysis of the Doran & Robbers model 
are only able to fit the distance to recombination well 
by changing f2 m and Hq from their true values. There- 
fore, including a strong, 1% prior on Hq in addition to 
the SN, CMB, and cluster data leads to nearly unbiased 
constraints on the curvature for this model of early dark 
energy. Including the full CMB constraints instead of just 
the distance priors would also help in this case, since the 
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2% early dark energy fraction at recombination would be 
detectable in Planck data. 

For models with massive neutrinos, the results of 
the growth reconstruction analysis are similar to those 
for early dark energy models. Massive neutrinos sup- 
press growth on scales below their free streaming length 
[6J, l65[ , affecting the growth evolution in a manner simi- 
lar to dark energy that tracks the dominant matter or 
radiation density at early times (u>oo = 0). Some of 
the strongest current cosmological limits on the neutrino 
mass come from the combination of the Lyman-a forest 
power spectrum with CMB data, with a 95% CL upper 
limit on the sum of neutrino masses of ^2 m v < 0.17 eV, 
assuming flat ACDM [66|. More conservative upper 
bounds limit the masses to J^ 171 ^ < 1-3 eV (95% CL) 
using the CMB alone and £ m v < 0.67 eV (95% CL) 
with CMB, BAO, and SN data (also assuming flat ACDM 
[l[ ) . On the other hand, neutrino oscillation experiments 
indicate that there must be at least one neutrino mass 
eigenstate with m u > 0.05 eV (e.g., [H [H). 

Using Woo as an effective parameter to absorb the ef- 
fects of massive neutrinos on distances and growth, the 
resulting curvature estimates tend to be biased toward 
high values of fix- The bias is similar to what we found 
for the Doran & Robbers early dark energy model above, 
and as in that case an additional Hq prior can reduce the 
curvature bias. Without any extra priors, the curvature 
can be biased by as much as ~ 0.5er even for the mini- 
mal mass of m v ~ 0.05 eV, so proper modeling of the 
effects of massive neutrinos is likely to be important for 
future curvature constraints from distances and growth. 

Using the sum of neutrino masses as an additional pa- 
rameter in the growth reconstruction analysis, instead of 
modeling the effects of neutrinos with Woo , produces un- 
biased curvature constraints. However, for models with 
massive neutrinos the uncertainty in curvature is fairly 
large due to a strong degeneracy between fix and ^ m v . 
In such a scenario, independent measurements of the neu- 
trino masses would greatly reduce the uncertainty in cur- 
vature from distance and growth probes. Future cosmo- 
logical measurements such as weak lensing of the CMB 
may be able to determine neutrino masses with an accu- 
racy of m v) ~ 0.05 eV, which should yield strong 
constraints on neutrino masses when combined with the 
results of terrestrial experiments (e.g., see [6^] for a re- 
view) . 

For the most general treatment of the high-redshift 
universe, one should ideally include parameters for both 
massive neutrinos and early dark energy models. We 
leave further study of this approach for future work, but 
note here that a simultaneous analysis of the impact of 
massive neutrinos and early dark energy on curvature 
estimation may be complicated by degeneracies between 
the two in the distance and growth observables. 



C. Dependence on SN and CMB data modeling 

Using the growth reconstruction method, we can study 
the relative contributions of scatter in the SN and CMB 
data to uncertainties in the reconstructed growth func- 
tion by only including the scatter in one of these data 
sets in the Monte Carlo simulations. For the predicted 
growth observables Pd(g|^K), the SN errors assumed for 
SNAP are the dominant source of uncertainty. More pre- 
cise CMB data than anticipated from Planck would have 
little impact on the uncertainties in the growth recon- 
struction. However, for the curvature constraint, more 
precise SN or CMB data would not significantly reduce 
<r(f2x) since the assumed 1% growth uncertainties for 
IXO clusters dominate the curvature uncertainty. 

If we could obtain more precise growth measurements 
with uncertainties less than a few tenths of a percent, 
then the curvature estimate would be limited by the pre- 
cision of the data from SNAP and Planck. However, the 
relative importance of these two data sets for curvature 
is opposite that for Pd(g|fix)- That is, given extremely 
precise growth measurements, cr(f2x) would be reduced 
more by improving CMB data beyond Planck than by re- 
ducing SN uncertainties. The reason for this is related to 
the covariance between Gq{z) and G(z max ): as explained 
above (see Fig. CMB data induce a positive correla- 
tion between the reconstructions of Gq(z) and G(z max ) 
from distance data. Since changing the curvature has op- 
posite effects on Gq(z) and G(z max ), the covariance from 
CMB constraints leads to more precise limits on curva- 
ture. Improving CMB constraints therefore reduces cur- 
vature uncertainty by strengthening the Gq(z) — G(z max ) 
correlation. 

The curvature forecasts are not strongly dependent on 
the redshift distribution of the supernova sample. For 
example, if we assume a uniform distribution of 2300 SNc 
over < z < 1.7 instead of the SNAP distribution of 
2000 SNe plus 300 low-z SNe, the curvature estimate 
from SNAP+Planck+IXO is unchanged. 

Interpolation of the distance-redshift relation between 
redshift bins or across gaps in the SN data introduces 
additional assumptions about the evolution of the ex- 
pansion rate. To test the impact of these assumptions, 
we have tried several distinct methods for interpolating 
between bins and found that the reconstructed growth 
changes by < 1%, and the mean value of fix in forecasts 
changes by < 5 x 10~ 4 . However, there is always an inher- 
ent assumption about smoothness of x( z ) m the growth 
reconstruction method that should be kept in mind when 
interpreting the results. 

D. Dependence on cluster data modeling 

We have assumed fairly optimistic future growth con- 
straints for the curvature forecasts, with 1% measure- 
ments of the growth function at z < 1.5 as anticipated 
for the proposed IXO cluster sample (HTJ] . We can eas- 



18 



ily study how the curvature forecasts depend on the as- 
sumed precision of the growth measurements by rescal- 
ing the Fisher matrix of the growth observablcs F g used 
in Eq. (f3~Tj) . We find that for less optimistic assump- 
tions about the growth uncertainties, it is still possi- 
ble to obtain interesting model-independent estimates of 
the curvature. For example, doubling the growth uncer- 
tainties to 2% increases the curvature uncertainty from 
a(n K ) = 0.0022 to a(n K ) = 0.0033, and for 3% growth 
measurements, (t(£Ik) = 0.0045. 

The forecasts presented in the previous sections as- 
sume that the cluster likelihood can be approximated 
as a multivariate Gaussian distribution for uncorrelated 
growth observables {Go(zi)} in Az = 0.1 redshift bins 
and G(z max ). In reality, we can expect that at least 
some of the growth observables will be correlated. Per- 
haps more importantly, the observed cluster abundance 
depends on distances and the expansion rate as well as 
the growth function, and there may be degeneracies be- 
tween these functions in the cluster likelihood. We might 
expect that the effect on curvature estimates of treating 
cluster abundance as purely a probe of the growth his- 
tory is small since the distance-redshift relation is well 
constrained by the SN and CMB data. Although a de- 
tailed study of curvature constraints using the full cluster 
likelihood is beyond the scope of this work, we present a 
simple test of this expectation here. 

To check how well the fiducial distance data sets con- 
strain degeneracies between distances and growth in the 
cluster likelihood, we compute the Fisher matrix for 
growth and distance assuming Poisson-distributed clus- 
ters [nUnl, 
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where 9 = ({Go(zj)}, {HoD(zj)}), taking derivatives of 
the number of clusters in each bin A(M,;, Zj) at the fidu- 
cial flat ACDM model. The mass bins {M{\ have width 
A InM ss 0.35 and start at a mass threshold M m - m , and 
the redshift bins {zj} have width Az = 0.1 and cover the 
range 0.1 < z < 1.5. 

As described in Sec. Ill C[ the abundance of massive 
clusters is exponentially sensitive to growth through the 
mass function dn/dM [Eq. ©]. However, the observed 
number of clusters depends on the cosmological model in 
other ways as well. In particular, the comoving volume 
element [Eq. ([8])] introduces additional dependence on 
the distance-redshift relation and expansion rate. The 
mass dependence of the effective volume in which a given 
survey probes clusters of a certain mass is also cosmology- 
dependent (and redshift-dependent) in general, although 
here for simplicity we take the volume to be constant 
above a mass threshold M min = 10 14 h~ x M Q . 

The main mass proxy proposed for IXO is the prod- 
uct of the X-ray temperature Tx and gas mass M gas , 
Yx = TA'-Mgas, which is expected to have a relati vely 
small, < 10% scatter based on cluster simulations U2 



Assuming lognormal scatter in the Yx — M relation, the 
total number of clusters in redshift bin Zj with width Az 
and mass bin Mi < M < M,_i_i is I73| 
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We include a 3% systematic error in M in uncorrelated 
Az = 0.1 redshift bins to represent the IXO forecast for 
the error in the normalization of the Yx — M relation 
from weak lensing mass measurements of 100 clusters per 
redshift bin [z|- 

Cluster mass estimates generally also depend on the as- 
sumed cosmology, although the exact dependence varies 
for different mass proxies. Recent results from simula- 
tions indicate that this additional cosmological depen- 
dence is fairly weak [74| , but it is nevertheless important 
to include in a full analysis of cluster data. For simplic- 
ity, however, we neglect the cosmological dependence of 
the cluster masses here and assume that the sensitivity 
of cluster abundances to cosmology is dominated by the 
mass function and the volume element. 

For the halo mass function [Eq. (Ol we use the 
paramctrization fit to simulations in Ref. [75j , 
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where A, a, and b are weakly redshift-dependent. As the 
authors of that paper note, it is probably more appro- 
priate for the wide variety of cosmological models con- 
sidered here to replace the redshift dependence of these 
parameters with dependence on the growth function. In 
fact, it may be important to include dependence on not 
only the instantaneous value of the growth function at 
the redshift of a cluster, but also the evolution of the 
growth function prior to that redshift [76$ . Recent stud- 
ies using N-body simulations or semi-analytic modeling 
of nonlinear growth in cosmologies with more complex 
dark energy evolution than ACDM have beg un to exam- 
ine such issues @ [ill HO, HI H3, HI, HI HE HI ■ For 
the approximate modeling of the cluster mass function 
here, however, we use the simpler redshift dependence 
given by Ref. [75l |. 

By taking the submatrix of F a p from Eq. (fM)) that 
corresponds to {Go(zj)} only and inverting to get the 
covariance matrix, we find that the growth uncertain- 
ties for fixed distances are ~ 1%. Marginalizing over the 
distances by inverting the full Fisher matrix, including 
both growth and distance variables, increases the growth 
function uncertainties by a factor of 5 — 10. However, 
adding SNAP-likc SN distance constraints with 1% accu- 
racy in {H D(zi)} reduces the growth uncertainties after 
marginalizing over distances to 1 — 3%. Recall from the 
beginning of this section that 3% growth function uncer- 
tainties approximately double the error in the curvature 
estimate relative to 1% growth uncertainties. 
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This test shows that distances at low z should be 
constrained well enough by the future SN and CMB 
data that the impact of additional cluster data can be 
roughly approximated by its expected accuracy on the 
growth function alone. Including degeneracies between 
growth and distances in the halo mass function and co- 
moving volume clement results in somewhat weaker but 
still interesting constraints on curvature. However, given 
the simplifying assumptions made here, curvature con- 
straints from distance and growth data using the full 
cluster likelihood are a subject that deserves further ex- 
ploration. 



E. Comparison with other methods 

The model-independent curvature constraint from dis- 
tances and growth has a forecasted accuracy of (t^k) = 
0.002 for SNAP SN, Planck CMB, and IXO cluster data, 
assuming that the true cosmology is close to fiat ACDM. 
This is very close to the accuracy expected in the model- 
dependent context of ACDM using the SN and CMB dis- 
tance data only. Thus the inclusion of growth informa- 
tion provides a curvature measurement that is free from 
possible biases due to assuming an incorrect form of the 
dark energy evolution without sacrificing precision. 

As we have seen, other types of measurements can play 
a similar role to the growth constraints. A 1% measure- 
ment of the Hubble constant and a 1% upper limit on the 
fraction of dark energy at recombination, when combined 
with SNAP SN and Planck CMB data, provide a model- 
independent curvature constraint with ~ 0.005. 
However, this constraint is strongly skewed with a long 
tail toward open models, and without growth informa- 
tion the constraints such data can place on cosmologies 
that have significant amounts of early dark energy are 
severely limited. 

The method of curvature estimation studied here 
can provide complementary constraints to other model- 
independent techniques. Knox [l!| proposed using pre- 
cise distance measurements, for example from BAO at 
high redshifts (z > 3), in comparison with the distance 
to recombination measured in CMB data to probe cur- 
vature. With future, percent-level BAO distances and 
Planck CMB data, this method is ex pec ted to attain an 
accuracy of a(Q K ) ~ 0.001 - 0.002 [lj, S3, S3- How- 
ever, this measurement depends on the assumption that 
the universe is matter-dominated between the rcdshift 
of the BAO measurement and recombination. As a re- 
sult, the estimated curvature is independent of the low- 
redshift dark energy modeling but still depends on the 
high-redshift dark energy evolution [l^|. This depen- 
dence is similar to the degeneracy between curvature and 
early dark energy that we find when comparing high-z 
SN distances to the CMB distance [Eq. (fl8])] . Additional 
information, such as a measurement of G(z) (or <j&), or 
BAO in the line-of-sight direction to probe H(z) at z ~ 3, 
could help to reduce the high-z model dependence of this 



method [lj|. 

The curvature measurement proposed by Bernstein 
fl8| is perhaps the most model-independent method since 
it relies only on the FRW form of the metric without 
assuming particular dark energy properties and is valid 
for alternative theories of gravity as well. Weak lensing 
galaxy-shear correlations can measure distance "trian- 
gles" involving the lens distance, source distance, and the 
distance between the two, and the relations between these 
distances are sensitive to curvature. Future lensing and 
galaxy surveys can provide a purely geometric test of cur- 
vature with an expected accuracy of <j{VLk) ~ 0.02 — 0.04 
[l8l [sol ] . This technique is likely to be even less model- 
dependent than the distance plus growth method de- 
scribed here, but the forecasted uncertainties are an order 
of magnitude larger. 

Given that each of these methods for obtaining model- 
independent curvature estimates relies on different types 
of data, it is difficult to compare forecasts directly. On 
the other hand, the existence of multiple methods using 
independent data sets means that there will be many op- 
portunities for cross-checks of the curvature estimates. 
In this sense, it is useful to have an array of model- 
independent methods available that will have different 
systematics, both theoretical and instrumental. 



VII. DISCUSSION 

With future supernova and CMB data sets making def- 
inite predictions for the growth of linear perturbations 
based on measured distances, precise measurements of 
the actual growth history will provide model-independent 
estimates of spatial curvature. Such constraints have im- 
portant implications for testing models of inflation, which 
predict J7k ~ 0, and for obtaining constraints on dark en- 
ergy models that are robust to uncertainty in curvature. 

If the true cosmology is similar to flat ACDM, as cur- 
rent data would suggest, then a combination of SNe 
from SNAP, CMB data from Planck, and X-ray clusters 
from IXO can measure the curvature with an accuracy 
of ct(Qk) ~ 0.002, making minimal assumptions about 
the dark energy evolution. The main constraint on cur- 
vature in this scenario comes from combining distance 
data with measurements of the evolution of the growth 
function at low redshifts. However, information about 
the normalization of growth relative to early times, ob- 
tained by comparing cluster abundances with the ampli- 
tude of CMB anisotropies, can significantly reduce the 
uncertainty in curvature. The extra information mainly 
comes from the covariance between the growth normal- 
ization and the low-rcdshift growth evolution required by 
the CMB constraint on the distance to recombination. 

These forecasts are robust to changes in the true value 
of the curvature and in the dark energy evolution at low 
redshift. However, if the true cosmology is significantly 
different from the concordance flat ACDM model at high 
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rcdshifts (z > 2). then the forecasted errors on fix in- 
crease by a factor of a few and the constraints depend 
more strongly on the normalization of the low-redshift 
growth function relative to early times. Although model- 
independent limits on curvature alone are weaker in this 
scenario, the combination of future SN, CMB, and cluster 
data would reduce the allowed model space to a narrow 
degeneracy between curvature and early dark energy or 
massive neutrinos. A precise independent measurement 
of the Hubble constant would mitigate possible biases in 
these constraints. 

Estimates of curvature from distances and growth are 
complementary to other model-independent techniques. 
The projected accuracy is similar to that expected from 
comparing distances measured at z ~ 3 to the distance to 
recombination [l9| , but the inclusion of growth informa- 
tion reduces dependence on assumptions about the high- 
redshift universe. Compared with proposed metric tests 
of curvature using weak lensing distance triangles [l8l ]. 
constraints from distances and growth are more model- 
dependent, in particular relying on the assumption that 
GR is valid on large scales and that dark energy cluster- 
ing does not significantly affect the growth observables. 
However, limits on curvature from combinations of dis- 
tance and growth data are potentially much more precise. 

We have focused here on the use of X-ray clusters to 
probe the growth history, but there are a number of 
other possible methods for measuring growth that may 
also provide interesting curvature constraints when com- 
bined with distance measurements. Correlations involv- 
ing weak lensing shear and galaxy density measurements 
should constrain both distances and growth in several 
redshift bins (e.g., 89]). Redshift space distortions of 
galaxy correlation functions can be used to constrain the 
growth rate rather than the integrated growth, so such 
data sets may provide interesting complementary con- 
straints to those from clusters or weak lensing. Fore- 
casts of the ability of these alternative probes of growth 
to measure curvature when combined with distance data 
are an interesting subject for future study. 

A fundamental assumption for using distances and 
growth to measure spatial curvature is that the relation 
between the growth history and the expansion rate is 
governed by GR. In the context of modified theories of 
gravity, many authors have studied the use of distances 
and growth to test for deviations from GR. These investi- 
gations typically make some simplifying assumptions re- 
garding spatial flatness and/or the dark energy evolution, 
and relaxing these assumptions could make such tests 
of gravity considerably more complicated. For example, 
at any particular scale the deviations in the distance- 
growth relation caused by modifying gravity may be dif- 
ficult to distinguish from the effects of nonzero curvature 
or dynamical dark energy. However, scale dependence of 
growth could still be a robust signature of certain classes 
of modified gravity theories even in the context of more 
general cosmologies. Regardless of our ability to use fu- 
ture data to distinguish between nonzero curvature, dy- 



namical dark energy, and modified gravity, finding hints 
of any of these possibilities would be an intriguing sign 
of physics beyond the standard cosmological model. 



Acknowledgments: The author would like to thank 
Wayne Hu, Cora Dvorkin, and Amol Upadhye for use- 
ful discussions. This work was supported by the NSF 
GRFP. 



APPENDIX A: GROWTH FUNCTION 
APPROXIMATION AT HIGH REDSHIFT 

At high redshift, assuming that the contributions of 
dark energy and curvature to H(z) are much smaller than 
that of the matter density, we can derive an approximate 
form for the growth history. In an Einstein-de Sitter uni- 
verse with matter only, the growth function is constant. 
We therefore start with an ansatz that the growth history 
in a universe with small fractions of dark energy and cur- 
vature is a small perturbation to the constant Einstein-de 
Sitter solution, expanding around some redshift Zi where 
G(zi) = Gj, 



G(z) = Gi + Cede 



C- 




(Al) 



where Cede <C Gj, C curv <C G^, and we have assumed 
that dark energy at early times can be parametrized by 
a constant effective equation of state Woo. The redshift 
dependence of the terms in Eq. (|A1[) is the same as that 
of Que{z) and flK(z) during matter domination. The 
growth solution must have this dependence in order to 
solve Eq. ^ since d In H/d In a contains terms propor- 
tional to fl-DE(z) and Qk(z). 

For this approximate solution to the growth equation, 
we neglect terms of second or higher order in Cede, 
C curv , ^de(^i)) an d ^k(zi)- In this limit, the quantity 
d In H/d In a appearing in Eq. ((6]) is 

dlxiH 3 1 — . , 3 _ . . , . „. 

— « -- + -n K (z) - -v^thmiz). (A2) 

Using Eqs. (|A1[) and (|A2[) in the differential equation 
for G(z) produces algebraic relations for the early dark 
energy and curvature coefficients: 



Cede 
C curv 



1 - W n 



Woo (5 - 6Woo) 

4 

-fl K {z l )G i . 



£luE(zi)Gi, 



7' 



(A3) 
(A4) 
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This approximation can be used to set the value of 
5'(z m&:x ) needed for the analytic growth reconstruction 
in Sec. R7J 



d\nG 



S(z B 



H ( £ max ) 



(1 + Z max ) Hq 



d\nG 



din a 



(A5) 



din a 



3(1 -Woo) , s 4 

•— "DE(,Z m axj — -"K(2maxJ- 

5 - 6w„o 7 



The approximation breaks down when either the cur- 
vature or early dark energy fraction is large. For 
^DE(zmax) < 0.25 and |Qk| < 0.01, the error in the 
approximation for S'/6 at z max ~ 1.5 is < 1%. The ap- 
proximate form of the growth rate in Eq. (|A5[) remains 
accurate even as Woo — ► despite the fact that the inte- 
grated growth function [Eq. (|A1|) ] becomes inaccurate at 
Woo > -1. 



APPENDIX B: MONTE CARLO SIMULATIONS 
OF GROWTH RECONSTRUCTION 

The procedure we use to estimate the covariance of 
growth reconstructed from distances (Sec. [Vj for future 
SN and CMB data is as follows: 

1. Define rcdshift bins {z^}, i = 1, ■■■,n z with zy = and 



Choose the assumed values of fix and early dark en- 
ergy/massive neutrino parameters (here, we use w^, as 
an example). 

Draw a realization of the SN data, specified by the red- 
shift distribution of Type la SNe and their magnitude 
errors: 



1 



Xa 



-s 



K 



) , (Bl) 



£stat,a ~t~ ^sys (Zq<)) 



where the statistical error e s tat,a is drawn for each 
SN (labeled by a) from a Gaussian with width 
CH D,a = 0.15-ffo-D(z Q )/(51oge), and the systematic 
error e sys (z Q ) is drawn for each Az = 0.1 redshift bin 
from a Gaussian with width 0.02[(1 + z)/2.7]/(5 loge). 
The SN distance-redshift relation is assumed to be un- 
biased relative to the fiducial model for the data, so 

(Ca> = 0. 

4. Estimate z(x) from the SN data by first estimating 
x(z) and then inverting the relation. We compute x( z ) 
in redshift bins Zs with width Az as 



X(zi) 



J2 a Xa exp[-(z Q 



2 /(2Az 2 



E a exp[-(z Q -z l ) 2 /(2Az 2 )] 



(B2) 



where z a and Xa are the redshifts and x values of SNe 
from step 3. The default bin width used here is Az — 
0.05. 

To invert this to obtain z(x), the estimated xi z i) 
must be monotonic. For large enough Az this is typi- 
cally not a problem since the fiducial x( z ) relation is al- 
ways monotonicaUy increasing, but even with wide red- 
shift bins there is a chance of having a few realizations 
with non-monotonic x( z i) estimates. We correct for 
this when it occurs by simply setting x( z i+i) = x( z i) 
in any bin for which Eq. (|B2j) gives x( z i+i) < x( z i)- 
The fact that the results are relatively independent of 
the redshift bin width indicates that this correction is 
not a significant source of error. 

SN coverage that is fairly uniform in z is important 
to avoid bias in the z{x) relation. For the anticipated 
SNAP distribution the largest biases are at the ends, 
z « and z ~ z max . The low- 2; bias has little im- 
pact on the growth reconstruction at z > 0.1, but 
the high-z bias can cause problems with connecting 
the reconstructed growth at z < z max to the fiducial 
z > Zmax growth history. To avoid such problems, we 
take z max for the growth reconstruction to be slightly 
smaller than the maximum SN redshift; in the case 
of SNAP where the SN distribution ends at z = 1.7, 
Zmax = 1-5 is sufficiently low to avoid problems relat- 
ing to bias in z(x). (Note that using a different value 
of z max results in slightly different definitions of early 
dark energy for the MCMC and growth reconstruction 
methods.) 

5. Estimate £ mM = dz / dx{z max ) using a linear fit to the 
high-z end of the z(x) relation. 

6. Draw simulated CMB distance data from the 2x2 
Fisher matrix for {ln^/Mpc), fi m /i 2 } (SecHTEj). Us- 
ing the assumed curvature and early dark energy pa- 
rameters in addition to the estimate of £ max from step 
5, set fi m and h to match the CMB constraints. 

7. Optionally include additional priors (e.g. Hq and dark 
energy fraction at last scattering). Each Monte Carlo 
simulation is weighted by the likelihood associated 
with these priors, and the entire set of simulations for a 
single {^k, Wqo} pair has a weight assigned to it equal 
to the average of the individual simulation weights. 

8. Compute the growth at z ma x relative to recombination, 
G(z max ), for the expansion history specified by J7k and 
Woo (step 2); E max (step 5); and f2 m and h (step 6). 

9. Given z(x) from step 4 and O m from step 5, iteratively 
solve the growth reconstruction equation, Eq. (|2T|) . For 
each iteration, 5' max is set to the average of its value 
in the previous iteration and the target value based on 
the approximate growth evolution at z > z max given 
by Eq. JM]). 
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10. Repeat steps 3 — 9 for many realizations of the SN and 
CMB data, and compute the mean and covariance (us- 
ing the weights from step 7) of the resulting estimates 
of g = {G(z max ), Go(zi)}. This produces the predicted 
growth observables at fixed fix and Woq. 



11. Repeat steps 2 — 10 for different curvature and 
early dark energy parameter values to compute 

Pd{g\^K,Woo). 
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